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ABSTRACT 



The long-term transitions of the black hole candidate Cygnus X-l (between the states 71, 
72, and 73) include the occasional appearance of a strong ~ MeV bump (71), whose strength 
appears to be anti-correlated with the continuum flux 400 keV) due to the Compton 
upscattering of cold disk photons by the inner, hot corona. We develop a self-consistent disk 
picture that accounts naturally for these transitions and their corresponding spectral variations. 
We argue that the bump is due to the self-Comptonization of bremsstrahlung photons emitted 
predominantly near the plane of the corona itself. Our results suggest that a decrease by a 
factor of ~ 2 in the viscosity parameter a is responsible for quenching this bump and driving 
the system to the 72 state, whereas a transition from 72 to 73 appears to be induced by an 
increase of about 25% in the accretion rate M. In view of the fact that most of the transitions 
observed in this source seem to be of the 72 — 73 variety, we conclude that much of the long 
term gamma-ray spectral variability in Cygnus X-l is due to these small fluctuations in M. The 
unusual appearance of the 71 state apparently reflects a change in the dissipative processes 
within the disk. 

Subject headings: accretion disks — black hole physics — radiation mechanisms: 
bremsstrahlung — radiation mechanisms: inverse Compton — relativity — stars: individual: 
Cygnus X-l 
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1. Introduction 

The black hole candidate Cygnus X-l has shown variability of its X-ray emission on 
all time scales ranging from milliseconds to years. Early observations of this strong hard 
X-ray source (e.g., Tananbaum et al. 1972; Holt et al. 1976; Ogawara et al. 1982) 
revealed long-term transitions between two principal states: the so-called "low-state" with 
a relatively small ratio of soft to hard X-ray flux, and the "high-state" corresponding to a 
softer spectrum. Generally speaking, the hard X-ray /7-ray emission has been successfully 
interpreted in terms of the Compton upscattering of soft disk photons by an optically thin, 
high-temperature plasma (~ 10 9 K) situated either at the inner edge of the disk (Thorne 
& Price 1975; Shapiro, Lightman & Eardley 1976), or within a coronal layer overlying 
the geometrically thin but optically thick cooler disk (Liang & Price 1977; Sunyaev & 
Titarchuk 1985). The first scenario could arise as a result of the secular instability present 
in the cool disk's inner region (Lightman & Eardley 1974), which would swell the optically 
thick, radiation-pressure dominated gas to a hot, gas-pressure dominated, optically thin 
cloud. 

More recent observations with the JPL High Resolution Gamma-Ray Spectroscopy ex- 
periment onboard HE AO 3 (Ling et al. 1987) have suggested a more complicated behavior 
than that assumed in a simplified two-state model. In particular, a new 7-ray emitting 
state of Cygnus X-l has been identified, in which the luminosity within an unusually strong 
7-ray "bump" between ~ 400 keV and ~ 2 MeV exceeds that of the 50 — 400 keV hard 
X-rays and is comparable to the overall emission below ~ 400 keV. This state, labeled 
71, together with the other low-state (dubbed 72) and the high-state (73), seem to form a 
sequence of transitions triggered, perhaps, by variations in the accretion rate, as proposed 
earlier for the two-state system (e.g., Kazanas 1986). 

Liang and Dermer (1988) have interpreted this 7-ray bump as due to the emission 
from a hot quasi-spherical pair-dominated cloud (distinct from the corona producing the 
continuum), in which the pair-balance condition determines the compactness parameter 
and the Thomson depth as a function of the equilibrium temperature. They find that 
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a reasonable fit to the "bump" may be obtained for a lepton temperature pa 400 keV 



for a black hole mass M. Our purpose in this Letter is to build on the strengths of 
these earlier models to arrive at a self-consistent accretion disk picture that accounts 
naturally for these transitions and their corresponding spectral variations. In particular, 
we relax the assumption that the inner hot plasma is uniform (Shapiro, Lightman & 
Eardley 1976; Liang & Dermer 1988) and consider a more realistic stratified structure 
consistent with gravitational settling. As we shall see, the secular instability still results 
in a two-temperature corona in the inner region of the disk, but the electron density in 
the disk plane may now be sufficiently high to produce a copious supply of bremsstrahlung 
photons. We shall argue that the self-Comptonization of this spectral component is the 
origin of the bump seen in the 71 state. Not surprisingly, the physical conditions in this 
inner hot corona are not unlike those inferred for the pair plasma by Liang and Dermer 
(1988), but here the lepton density is due mostly to accretion, not pair production. 



Our starting point is to assume a standard a-disk model (Shakura & Sunyaev 1973), 
which however is subject to the secular instability in the inner region that drives the 
disk from a cool state (with electron temperature T e ~ 10 6 K) to the hot, two-fluid 
state (T e ~ 10 9 K and ion temperature Tj 3> T e ; Shapiro, Lightman & Eardley 1976, 
hereafter SLE). Unlike the earlier work, we do not employ vertically-averaged equations 
for the coronal structure, but rather allow for density stratification as discussed above. For 
simplicity, we assume constant temperatures in the vertical (i.e., z) direction, though not 
in the (cylindrical) radial direction (R). As such, SLE's equations (3) and (4) are replaced 



and a pair cloud radius of about 9 r g , where r g = 2GM/c 2 is the Schwarzschild radius 



2. A Self-Consistent Two- Temperature Disk 



2.1. The Structure Equations 
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respectively, where r 2 = R 2 + z 2 , P g is the gas pressure, p is the density and J = 1 — 
{Zrg/r) 1 / 2 . We ignore the contribution to the lepton number density from pair production 
and check a posteriori that this is indeed a valid approximation. In addition, the electron- 
ion energy exchange flux (SLE equation 10) generalizes to 



J dF = J ^v E pk(T t - T e )/m p dz , 



(3) 



where ve is the electron-ion coupling rate. 

It is not difficult to see that with these modifications, the two-temperature structure 
of the hot inner corona is now specified by 



p = po exp(-z 2 /zQ) , 



(4) 
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and the temperatures are given by 
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This prescription, however, is inappropriate whenever the physical conditions are such 
that the protons attain sufficient energy to escape from the system. Since the virial tem- 
perature 

I - / r _ \ 

(8) 



T = -m c 2 (Tg 



where m p is the proton mass, decreases with increasing radius r, it is anticipated that the 
ion temperature Tj may exceed T vir at large radii. In that case, a coronal mass (and energy) 
outflow will prevent Tj from greatly surpassing T V i r . For the purpose of this calculation, 
we shall assume that the fraction of dissipated power lost in this wind is at most a small 
fraction of the radiated luminosity. (This restriction will be relaxed in future radiative- 
hydro dynamical simulations of the accretion geometry.) To simplify matters further, we 
shall also assume that in this region (where the upper layers of the corona merge into a 
"wind"), the ion temperature is in fact equal to the virial temperature (though of course 
this can only be approximately correct). Instead of the above set of equations, the relevant 
expressions that determine the coronal structure are then Equations (1), (3), and the 
condition Tj = T vir , which yield 

p = Po (i+. 2 AR 2 r 3/2 , (9) 

and 

Te = 5 - 2xl ° 9K (-Vfinrm^ 3 

\r g J \5 x 10 5 gem 6 
In writing these equations, we have assumed that although the outer coronal layers may 
merge into a transonic flow (at say z son i c ), most of the structure of interest to us here lies 
at z <C z sonic, for which Equation (1) is still an adequate representation of the vertical 
density profile near the equatorial plane. 

2.2. Method Of Solution 

Following SLE, we have taken the outer boundary of the two-temperature inner region 
to lie at a radius R given by the condition P r {R Q ) = 3P g (R ), where P r is the radiation 
pressure in the a-disk. We note, however, that our results are insensitive to the actual 
location of this boundary, as long as P r (Ro) ~ 0[P g (Ro)]- For example, when Rq is instead 
fixed by the condition P r (Ro) = 3/2 P g (Ro), our spectra differ from those presented below 
by at most 10%, well within the observational uncertainties. In this model, the cooling 
results from the inverse Comptonization of both the cold disk photons penetrating into the 
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corona and the bremsstrahlung radiation produced in this inner hot region. To handle this 
dichotomy, we divide the corona into concentric rings and proceed as follows. We assume 
a fractional cooling by each spectral component in the first outer zone and then iterate 
over the mid-plane density until the local bremsstrahlung-self-Compton emissivity (within 
this ring) is balanced by the assumed fraction of the dissipation rate at this location. 
Within any given ring, the contribution to the overall flux from the Comptonized (outer- 
disk) photons is proportional to the total number of scattering particles and the electron 
temperature within that ring, and the attenuated soft photon number density, as reflected 
by the dependence of the Kompane'ets equation on these parameters. Thus, in subsequent 
rings, we iterate on the local mid-plane density until the calculated fraction agrees with 
the value extrapolated on the basis of this proportionality from the adjacent outer ring. 
From this we determine the overall coronal luminosity, including both components, and 
then iterate on the assumed fraction in the outer ring until this power matches the total 
dissipation rate inside the corona. We check ring by ring whether the first set of structure 
equations (see §2.1. above) yields an ion temperature Tj > T v i r . If mass loss is indicated 
at a particular radius, we then use the second set of equations to determine the local 
structure. 

The cooling rate (and concomitantly the spectrum) is obtained by solving the relativis- 
ts steady-state Fokker-Planck equation, whose diffusion coefficient is valid for arbitrary 
photon and electron energies (Prasad et al. 1988). At the temperatures of interest in our 
problem, the use of the non-relativistic Kompane'ets equation (Kompane'ets 1956) is not 
valid, which in the past has led to the development of a (time-consuming) Monte-Carlo 
scheme in order to treat the inverse Compton scattering correctly (e.g., Liang & Dermer 
1988). With the recent development of the exact analytical formula for this (relativistically- 
correct) diffusion coefficient of the Compton Fokker-Planck equation, we are now able to 
circumvent this difficulty by simply solving the relativistic Kompane'ets equation without 
any approximation. In this regard, our approach differs from that of SLE in at least 2 ways. 
First, our source term includes both the cool disk photons penetrating the y & 1 region of 
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the corona and the bremsstrahlung photons produced within the corona itself (predomi- 
nantly near the plane). Here, y = (4/cT e /m e c 2 )Max(r, r 2 ) is the dimensionless parameter 
that characterizes Comptonization, and r is the scattering optical depth. Secondly, since 
our equation is valid at all energies, our results are not subject to the restrictions (e.g., 
in T e ) and errors inherent in Cooper's modification to the Kompane'ets equation (Cooper 
1971). Of course, the main result of our calculations is the spectrum, and this is deter- 
mined along the lines developed by SLE, in which the energy-dependent rate of photon 
loss appearing in the Fokker-Planck equation also acts as the source term for the energy- 
dependent photon flux. 

3. The Gamma- Ray Spectral Variability 

Since the structure of the hot inner corona is here calculated self-consistently, we have 
at most only 3 physical quantities that need to be fixed. These are the black hole mass M, 
the accretion rate M, and the viscosity parameter a. A very conservative lower mass limit 
for Cygnus X-l is 3.4 M & (e.g., Paczynski 1974), whereas the most likely value is thought 
to lie in the range 9 — 15 M (e.g., Avni & Bahcall 1975). We have chosen to work with 
a value M = 10 M , though again our results are insensitive to the actual mass as long 
as 5 M & M & 15 M . The basis for this "robustness" is the fact that distances scale as 
the Schwarzschild radius r g = 2GM/c 2 , so that the overall energy budget (~ GMM/r) is 
independent of M (to first order) for a given M. Nonetheless, the coronal structure does 
change with M, but we have found only ^ 10% differences in the spectra over this mass 
range and have therefore not considered M to be a principal parameter characterizing our 
solutions. In addition, although we expect (slight) variations in M between the various 
states of Cygnus X-l, the gross X-ray /7-ray energetics for an assumed source distance of 
2.5 kpc restrict the accretion rate to a value ps 10 18 g s -1 . 

Our solutions for the 71, 72, and 73 states are compared to the data in Figures 1, 2, 
and 3. The data used by SLE overlap most significantly with the 71 spectrum and have 
therefore been included in Figure 1. We note also that we have excluded the 1.5 MeV point 
in the HE AO 3 data since it may be contaminated by a background line (Liang and Dermer 
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1988). In each case, the dashed curve shows the contribution from the bremsstrahlung/self- 
Comptonized (BSC) photons originating from within the corona itself, whereas the solid 
curve represents the overall spectrum comprising the BSC radiation, the cold disk photons 
that are Compton upscattered in the y & 1 region of the corona, and the thermal emission 
from the disk itself (which dominates at E ^ 3 keV). What emerges from these results 
is the interesting correlation between the variation in a and the gradual shift of emitted 
power from the 7-ray continuum into the broad ~ MeV bump (see Figure 4). In particular, 
a factor 2 decrease in a is sufficient to drive the system from the 71 to the 72 states. Thus, 
the emergence of the bump in the 72 and especially the 71 states must be accompanied by 
a corresponding reduction in the continuum flux below about 400 keV, as confirmed by 
the data. On the other hand, an increase in M raises the flux in both the continuum and 
the bump. It appears that an increase of about 25% in M is responsible for driving the 
system from the 72 to the 73 states (see Figure 5). 

Figures 6, 7, 8, and 9 demonstrate how the internal structure of the corona changes 
in response to a variation in either a (characterizing the transition from the 71 to the 72 
states) or the accretion rate (which induces the transition from 72 to 73). In Figure 6, 
the most significant aspect of the electron temperature profile is that the corona becomes 
optically thick toward small radii due to gravitational stratification (see the profiles of 
mid-plane density shown in Figure 7), and T e increases rapidly in this region to enhance 
the bremsstrahlung-self-Compton cooling process. The profiles of ion temperature are 
shown in Figure 8. In all three states, Tj increases with radius until it reaches the virial 
temperature T vir , though the radius at which this occurs is significantly larger in the 71 
state than the other two. The impact of this change in the ion temperature is reflected 
in the profiles of the y-parameter shown in Figure 9. Here, the solid curves correspond to 
the location in the corona where y = 1, whereas the dashed curves are for y — 0.5. As one 
would expect, the corona appears more extended in those regions where Tj pa T vir . The 
temperature, T e jumps to a value ss 50 keV near the boundary (at Ro ~ 15r s ) where the 
instability first arises, and then continues to increase towards 350 keV as the matter 
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approaches ~ 6 — 7r g . Thus, since most of the bremsstrahlung photons are emitted at 
y > 1, it is evident from Figure 5 that the location of the bump at ~ 1 MeV is a natural 
consequence of the self-Comptonization of the radiation, which approaches a Wien peak 
centered at roughly 3kT e in this inner hot region. However, the midplane density po is 
lower in the 72 and 73 states compared to 71 , and since the bremsstrahlung intensity scales 
as Po, the appearance of the bump is more evident in 71 than 72 and 73. 

Using these results, we may now justify a posteriori our neglect of pairs in the equi- 
librium lepton number density. Under the simplifying (and conservative) assumption that 
the radiation distribution within the corona is isotropic, we estimate the pair production 
rate via 77 interactions to be about 9 x 10 40 s _1 in the 71 state. (The pair creation rate 
is smaller in the 72 and 73 states due to the relatively lower number density of E > 511 
keV 7-rays.) The average (accretion) electron number density is (n e ) « 5 x 10 17 cm -3 , 
for which the positron (in-flight) annihilation time is t a = l/(?i e cry c) «lx 10~ 4 s. This 
is in fact the shortest time scale associated with the positron distribution, providing an 
upper limit to the positron number density n pa i rs . We therefore expect the positron aver- 
age equilibrium number density to be (n pa i rs ) ^ 8 x 10 13 cm -3 for the coronal structure 
depicted in Figure 5, in which the scale size is 3 x 10 7 cm. As such, n pa i rs <C n e and 
our approach would appear to be valid, at least in the case of Cygnus X-l. Incidentally, if 
we further assume that the subsequent conversion of these pairs into 511 keV photons is 
an efficient one, the 511 keV luminosity should be ss 2.8 x 10 36 ergs s -1 , or about 10 times 
weaker than the luminosity in the bump. 

This work was supported in part by NSF grant PHY 88-57218, NASA grant NAGW- 
2380, and the Alfred P. Sloan Foundation. We are grateful to Edison Liang for a helpful 
suggestion, and to the anonymous referee for his thoughtful comments and for pointing 
out an error in the original presentation of our spectra. 
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Figure Captions 

Figure 1. Theoretical spectrum corresponding to the 71 state of Cygnus X-l. Dashed 
curve: bremsstrahlung/self-Compton component (BSC); Solid curve: overall spectrum, 
including the BSC component, the disk radiation upscattered within the y & 1 region of 
the inner corona, and the thermal disk radiation at ^ 3 keV. Data for E <S 50 keV are 
taken from Ling et al. 1987. Data below ~ 50 keV (which however were collected at a 
different epoc) are from Schreier et al. 1971. 

Figure 2. Same as Figure 1, except now for the 72 state of Cygnus X-l. All data are 
from Ling et al. 1987. 

Figure 3. Same as Figure 1, except now for the 73 state of Cygnus X-l. All data are 
from Ling et al. 1987. 

Figure 4. Comparison of the 71 and 72 theoretical spectra, showing the cross-over at 
approximately 400 keV, and the anti-correlation between the strength of the bump at as 1 
MeV and the continuum flux below this cross-over energy. 

Figure 5. Comparison of the 72 and 73 theoretical spectra, showing the overall shift in 
flux, for both the continuum and the bump, when M changes by ~ 25%. 

Figure 6. The run of electron temperature T e as a function of the radius r in the disk. 
The onset of the secular instability at r a* 15 — 17 r g is signaled by a rapid increase in 
both T e (to <S 50 keV), as specified by Equation (6) or (10). Because of gravitational 
stratification, the corona becomes optically thick toward smaller radii, and T e increases 
further to enhance the cooling rate due to bremsstrahlung-self-Compton emission in this 
region. 



Figure 7. Profiles of the mid-plane density as a function of radius for each of the three 
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gamma-ray states. The slight jump in density reflects a transition from the bound corona 
(for which Tj < T vir ) to a "coronal wind" zone (in which Tj « T vir ). 

Figure 8. The run of ion temperature T; as a function of radius for each of the three 
gamma-ray states. Tj is limited from above by the value of the virial temperature T vir . 

Figure 9. Contours of the Compton y-parameter for y = 0.5 (doted) and 1.0 (solid), (a) 
71 state, (b) 72 state, and (c) 73 state. Distances are in units of the Schwarzschild radius 
r g = 2GM/c 2 . 
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